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Abstract 

We describe a practical method for constructing a nontrivial homomorphism be- 
tween two Verma modules of an arbitrary semisimple Lie algebra. With some addi- 
tions the method generalises to the affine case. 

A theorem of Verma, Bernstein-Gerfand-Gel'fand gives a straightforward criterion for 
the existence of a nontrivial homomorphism between Verma modules. Moreover, the the- 
orem states that such homomorphisms are always injective. In this paper we consider the 
problem of explicitly constructing such a homomorphism if it exists. This boils down to 
constructing a certain element in the universal enveloping algebra of the negative part of 
the semismiple Lie algebra. 

There are several methods known to solve this problem. Firstly, one can try and find 
explicit formulas. In this approach one fixes the type (but not the rank). This has been 
carried out for type An in Section 5, and for the similar problem in the quantum 
group case in H] , |5] , . In |1] root systems of all types are considered, and the solution is 
given relative to so-called straight roots, using a special basis of the universal enveloping 
algebra (not of Poincare-Birkhoff-Witt type). In [2], [HI the solution is given for types An 
and Dn for all roots, in a Poincare-Birkhoff-Witt basis. Our approach compares to this in 
the sense that we have an algorithm that, given any root of a fixed root system, computes 
a general formula relative to any given Poincare-Birkhoff-Witt basis (see Section E)). 

A second approach is described in ^Hj, which gives a general construction of homo- 
morphisms between Verma modules. However, it is not easy to see how to carry out this 
construction in practice. The method described here is a variant of the construction in 
|18j . the difference being that we are able to obtain the homomorphism explicitly. 

In Section^ of this paper we review the theoretical concepts and notation that we use, 
and describe the problem we deal with. In Section |21 we derive a few commutation formulas 
in the field of fractions of f/(n^). Then in Section |21 the construction of a homomorphism 
between Verma modules is described. In Section H] we briefly comment on the problem of 
finding compositions of inclusions. In Section El we comment on the analogous problem for 
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affine algebras, and we show how our algorithm generalises to that case. Finally in SectionlHl 
we give an application of the algorithm to the problem of constructing irreducible modules. 
This is based on a result by P. Littelmann. 

I have implemented the algorithms described in this paper in the computer algebra 
system GAP4 (jZI). Sections El and IHl contain tables of running times. All computations for 
these have been done on a PII 600 Mhz processor, with lOOM of memory of GAP. 

1 Preliminaries 

Let g be a semisimple Lie algebra, with root system $, relative to a Cartan subalgebra f). 
We let A = {ai, . . . , ai} be a fixed set of simple roots. Let = {ai, . . . , as} be the set of 
positive roots (note that here the simple roots are listed first). Then there are root vectors 
Ui = X-ai, Xi = (for I < i < s) , and basis vectors hi E I) (for I < i < I), such that the 
set {xi, . . . ,Xs,yi, . . . ,ys,hi, . . . , hi} forms a Chevalley basis of q (cf. [IHl)- We have that 
= © () © n"*", where n~, n"*" are the subalgebras spanned by the Ui, Xi respectively. 
In the sequel, if /5 = aj G then we also write in place of 

We let P denote the integral weight lattice spanned by the fundamental weights Ai, . . . , A;. 
Also QP = QAi + ■ ■ ■ + QA,. For A, G QP we write /i<Aif/i = A- ^'^^ kiOi, where 
ki G Z>o. Then < is a partial order on QP. 

For o; G $ we have the reflection '■ QP — ^ QP, given by Sq(A) = A — (A, 

Let U (g) denote the universal enveloping algebra of g. We consider U (g) as a g-module 
by left multiplication. Let A = ^ ajAj G QP, and let J(A) be the g-submodule of U (g) 
generated by /ij — + 1 for 1 < i < / and Xi for 1 < i < s. Then M(A) = f/(g)/J(A) is 
a g-module. It is called a Verma module. As f/(g) = U{n^) © J(A) we see that U{n^) = 
M(A) (as f/(n~)-modules). Let vx denote the image of 1 under this isomorphism. Then 
hi-vx = [tti — l)f A, and Xi-vx = 0. Furthermore, all other elements of M(A) can be written 
as Y ■ Vx, where Y G U{n~). 

Let u = Yl\=i hc^i, where G Z>q. Then we let f/(n^)y be the span of all yt-^ ■ ■ - yi^ 
such that + ■ ■ ■ + = u. 

For a proof of the following theorem we refer to P], [3]. 

Theorem 1 (Verma, Bernstein-Gel'fand-Gel'fand) Let A, G QP, and set 

P^,A = Hom^(g)(M(/i),M(A)). 

Then 

1. dimP^^A < 1, 

2. non-trivial elements of R^^x (^^e injective, 

3. dim P^^A = 1 ^/ one? only if there are positive roots , . . . , ctj^ such that 

< Sa,, (/i) < Sa^^Sa,^ (^) < " " " < S^,^ ■ ■ ■ S^,^ (/x) = A. 
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The problem we consider is to construct a non-trivial element in i?^ A if dimi?^ A = 1- 
By Theorem ^ this boils down to finding an element in i?^ ^ if /i = Sq(A) = A — (A, Q^)a 
and (A, a^) G Z>o. Suppose that we are in this situation, and set h = (A,a^). An 
element Y ■ vx & M{X), where Y G U{n~) is said to be singular if Xa ■ {Y ■ vx) = for 
a G Let G -R^,a be a non-trivial ?7(g)-homomorphism. Then ip{Vfj,) = Y ■ vx for 
some Y G with Y ■ vx singular. We have hiy = yhi — (z/, a^)?/ for all y G f/(n^)j,. 

Hence hi ■ {y ■ vx) = {{X — u, a^) — l)yvx. So, as hi ■ (Yvx) = ((/i, a^) — l)Yvx we see that 
Y G U{n~)ha- Conversely, if we have a F G U{n~)ha such that Y ■ vx is singular, then 
ip '■ M{fi) — ^ M(A) defined by = Y'Y -vx will be a non-trivial element of -R^,a- So 

the problem reduces to finding a F G U{y\r)ha such that F ■ f a is singular. Note that this 
can be done by writing down a basis for U{y\r)ha and computing a set of linear equations 
for Y . However, this algorithm becomes rather cumbersome if dimf/(n~)/ja gets large. We 
will describe a more direct method. 

2 The field of fractions 

From 3J, §3.6 we recall that U{vc) has a (non- commutative) field of fractions, denoted by 
K{n~). It consists of all elements ab^^ for a G f/(n~), h G U{x\r) \ {0}. For the definitions 
of addition and multiplication in K{x\~) we refer to l3J, §3.6. They imply aa~^ = a~^a = 1. 

Let a,P E If a + (3 G $^ then we let Na^p be the scalar such that [yayVis] = 

—Na^/3ya+i3- Also Set Pa^i3 = {ia + il3 I j > 0} n Then there are seven possibilities 



for 




(I) Po.,P 




(11) Pa,/3 




(ni) 


= {a, /5, a + /?, a -|- 2/5}, 


(IV) 


= {a, /5, a + /?, 2a + /?}, 


(V) 


= {a, /5, a + /3, 2a + /3, 3a + /3, 3a + 2/5} 


(VI) P.,p 


= {a, /5, a + /?, a + 2/?, a + 3/?, 2a + 3/3}, 


(VII) 


= {a, /5, a + /?, 2a + /?, a + 2/?}. 



Lemma 2 /n case (I) we have y^y^ = VaV^ for all m, n G Z. 

Proof. If n > then yisy^ = VaUi^- Multiplying this relation on the left and on the right 
by Va"' get ypy^^ = ya^yp- So we have ypy^ = UaUp fo^ all n G Z. From this it follows 
that y^y^ = y^y^ for m > 0, n G Z. If we now multiply this from the left and the right 
by y^'^ we get the result for m < as well. □ 



Since 
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these binomial coefficients are defined for arbitrary n G Q, and k G Z>o. In fact, we see 
that (2) is a polynomial of degree k in n. Note also that if n G Z and < n < k then the 
coefficient is 0. 

Lemma 3 In case (II) we have for m > 0, n G Z, 



Proof. First of all, this formula is known for m,n > (see, e.g., [H]). In particular, for 
n > we have ?//3?/2 = + Na,i3ny^~^ya+i3- If we multiply this relation on the left and 
the right with ?/~", and use Lemma El then we get it for all n G Z. Now the formula for 
m > 1 is proved by induction. □ 

Lemma 4 In case (III) we have for m > 0, n G Z, 

„.m n _ m,n n-k-l ■m-k-21 k I 

yp ya — 2-^ ^k,l ya yp ya+l3ya+2f3^ 

k,l>0 
k+2l<m 



where 

^k,i -^''^^P^2''P^-+i'> \k + l)\k + 2l)\ I 
Proof. This goes in exactly the same way as the proof of Lemma El Q 

Lemma 5 In case (IV) we have for m > 0, n & Z, 

..■m n _ m,n n~k-2l m-k~l k I 

yp ya — 2-^ ^k,l ya y/3 ya+(3y2a+f3^ 



k,l>0 
k+Km 



where 



Proof. Again we get the formula for m, n > from jUj. In this case the formula for m = 1, 
n > reads 

.n _ „.n„, I AT- „„.n-l„, , AT AT / \„.n-2„ 



ypya = y^yp + Na.pnyl y^+p + Nc,pN^^^+p y^j yl y2a+(3- 

If we multiply this on the left and the right by ?/~", and use Lemmas EUni we get the same 
relation with n replaced by —n. So the case m = 1, G Z follows. The formula for m > 1 
now follows by induction. □ 
The cases (V), (VI), (VII) can only occur when the root system has a component of 
type 6*2 • We omit the formulas for these cases; they can easily be derived from those given 
in 0. 
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Now let a = y^^ ■ ■ -y^" be a monomial in U{n ). For /3 G $^ and m,n E "Z consider 
the element y^ay^^. By repeatedly applying Lemmas El IH El "we see that 

y^ay^- = Yl • • • ' fc*)l/r""''''"'""'*'*«(^i' • • • ' ^*)- (1) 

(fcl,.--,fct)G/ 

Here the a(/ci, . . . ,kt) G t/(n^), the (finite) index set J, the Pi G Z>o are all independent 
of n, they only depend on a. Only the exponents of and the coefficients c{ki, . . . ,kt) 
(which are polynomials in n) depend on n. 

Now we take m, n G Q such that m — n E Z. Then we define y^ay^^ to be the right- 
hand side of ((H), and we say that y'^'ay'^'^ is an element of K{n^). More generally, if Y is 
a linear combination of monomials, and m, n G Q such that m — n E Z then y^^Yy^^ is an 
element of K{n~). 

3 Constructing singular vectors 

Here we suppose that we are given a A G QP and a G $^ with (A, a^) = /i G Z>o- The 
problem is to find a F G U{n^)ha such that F ■ f a is a singular vector. 

We recall that / = |A| is the rank of the root system. Let 1 < i < I, then 

Xiyl ■ vx = r((A, a^) - r)y\'^ ■ vx- (2) 

Lemma 6 Suppose that a G A, i.e., a = ai, I < i < I. Then y^ ■ vx is a singular vector. 

Proof. This follows from ((21), cf. the proof of Lemma 2. □ 
Note that this solves the problem when q = sl2- So in the remainder we will assume 
that the rank of the root system is at least 2. By an embedding : M{fi) ^ M (A) we will 
always mean an injective ?7(g)-homomorphism. 

Lemma 7 Suppose that z/, G P, and f3 E A is such that m = (z/, is a non-negative 
integer. Suppose further that we have an embedding ip : M(z/) M{r]) given by ip{v^) = 
Yvrf. Set n = {rj,(3^). Then yJ^Yy^"' is an element ofU{n~) and we have an embedding 
: M{spiy) ^ M{spr]) given by 0(fs^,.) = V^Yy^"" ■ Vs^r,- 

Proof. If n < then the first statement is clear. The embedding is the composition 
Mlspu) "—>■ M{v) ^ M{r]) ^ Mlsprj), where the first and the third maps follow from 
Lemma (HI 

If n > 0, then we view M^sprj) as a submodule of Mijj). We have f^^^ = y^Vr, 
(Lemma ini). Set v = y^Yvn', then f is a singular vector (being the image of Vg/^u under 
M(s/3i/) ^ M(i/) ^ Mijj)). We claim that v G M^sprj). Suppose that this claim is proved. 
Then there is a F' G U{n') such that v = Y'vspn- ^^^^ means that y'^Y = Y'y^, and 
the lemma follows. 

The claim above is proved in T^. For the sake of completeness we transcribe the 
argument. Set V = Mirf) / M^sprf), and let Vy denote the image of 'ipivy) in V] then 
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V,, = X ■ v,j, for some X G U{n~). For A; > write y'^X = Xiy^\ By increasing k we 
can get ki arbitrarily large (cf. 0, Lemma 7.6.9; it also follows by straightforward weight 
considerations). By Lemma El we know that y^Vj^ G M^Sjsrj). Therefore there is a A; > 
such that y^Vy = 0. Then by using Q we see that the smallest such k must be equal to 
m. □ 

Proposition 8 Let z/, ?7 G QP be such that v = s^{ri) = rj — k'-f, where 7 G and k G Z>o- 
Let Y G f/(n^)fc-y be such that Y ■ Vr, is singular. Let /? G A, /5 7^ 7. Set m = (z/, 
n = {ri,(3'^). Then y'^Yy^^ is an element of K{n~); it is even an element of U{n~). 
Secondly, we have an embedding (p : Mi^spv) ^ M^sprj) given by (p^Vs/^u) = y^'^Vp^ " '^sprj- 

Proof. We have that m — n = —k{'y,P'^) G Z, so y^'Yy^"' is an element of K{n~). 

Set = {/i G QP I (/i, 7^) = k}, which is a hyperplane in QP, containing 77. Let 
{oi, . . . , at} be a basis of U{n~)k-y Take /i = Yl\=i "^i^i ^ ^ ^"^^ set /i = = — A;7. 

Then by Theorem Q there is a = X]i=i Ci'^i such that Y^ ■ is singular. Here the Q are 
polynomial functions of the rj. (Indeed, if we write linear equations for the Y^, then the 
coefficients depend linearly on the r^. Hence the coefficients of a solution are polynomial 
functions of the r^.) 

Set p = (/i, q = (/i, Then Y' = y^Yfj_y^'' = J2j (^j^j^ where the bj are linearly 
independent elements of K{n~), and the cj are coefficients that depend polynomially on 
the Tj. Now Lemma [7| implies that if the G Z and p > 0, then Y' G U{n~). Let now 
j be such that bj ^ U{n~). If the r-j G Z and p > 0, then Cj = 0. Suppose that /? = ajj,, 
the ZQ-th simple root. Then the condition p > amounts to > k{'j,(3'^). We have that 
fi E V if and only if Yl\=i "^i^i ~ where the Ui are certain elements of Z. Also, since 
/3 7^ 7 at least one ui 7^ with i 7^ iq. We see that the requirement > /c(7,/5^) cuts a 
half space FT off V . Furthermore fl P is an (/ — 1) -dimensional lattice in V (cf. [1 ). The 
conclusion is that Cj = if G fl P. Since the Cj are polynomials in the r^, it follows 
that Cj = if /i G V. In particular, y^Yy^"" lies in [/(n~). 

Finally we note that Y'-Vg^^^ is singular, by the same arguments. (Indeed, Xi-(Y' -Vsg^) = 
J2j fj^j ■ "^s^M where the fj are polynomials in the rj, and the Zj are elements of U{n^). 
Since the are zero when /i G fl P we have that fj = when /i G V^.) In particular, 
yfYy^"^ ■ Vs^n is singular. □ 

Example 9 To illustrate the argument in the preceding proof, consider the Lie algebra 
of type As, with simple roots a,/?, 7 (with [3 corresponding to the middle node of the 
Dynkin diagram). Then it is possible to choose a Chevalley basis such that [?/a, yp] = ya+is, 

[ya,yf3+^] = ya+ci+j, b/3,1/7] = [yi,ya+f}] = -ya+p+^- Set ai = ?/«i//3i/7, = y-^ya+p, 

c^s = yayp+'j, Ci4 = |/«+/3+7- Then {oi, 02, as, 04} is a basis of U{n~)a+i3+'y 

We abbreviate a weight riAi + r2\2 + r-^Xs by (ri,r2,rs). Let be the hyperplane in 

QP consisting of all weights /x such that (/x, {a + P + 7)^) = 1, i.e., V = {{ri,r2,rs) \ 

ri + r2 + rs = 1}. Let /i = (ri, r2, r^) G K and set /i = Sa+/3+7(/i) = (n - 1, r2, rg - 1). Set 
= «! — ria2 — (ri + r2)a3 — rir^a^; then ■ is singular. Set p = {ft, a^) = ri — 1 and 

q= =ri. Now F' = ?/Py^?/-<? = y^y^ - (^ +r2)y/3+^ + ri(l -ri -r2 -rs)^/^^l/c.+/^+7• 
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According to Lemma [7| this is an element of U{n^) whenever (ri,r2,r3) G V with the 
integral and p > 0. Therefore the coefficient of y^^ya+/3+'y has to vanish, which is indeed 
the case. We see that Y' lies in U{n~) for all (ri, r2,r3) e V. 

Now we return to the situation of the beginning of the section. We have A G QP, 
a G with (A,a;^) = h E Z>o- Set yU = Sa(A) = A — /la. To obtain an embedding 
M{fi) M{X), we perform the following steps: 

1. Select . . . , /5r G A and positive roots ao, . . . ,ar in the following way. Set = a, 
and k = 0. Then: 

(a) If ak G A, then set r = A; and go to step 2. 

(b) Otherwise, let Pk+i G A be such that {ak,P^_^_i) > 0, and set ak+i = 5/3^.^-^(0^), 
and k := k + 1. Return to (a). 

2. Set /3 = a,. G A. For 1 < k < r set ak = —{ii,sp^---Sf3^_^{[3ky), and hk = 

3. Set Yo = y% and for < A; < r - 1: 

n+i = ylllYkvlzl. 

Remark. Note that the (3k+i in step 1 (b) exists because otherwise (afc,7^) < for all 
7 G A, and this implies that the set A U {cxk} is linearly independent (cf. jlHj, Chapter 
IV, Lemma 1), which is not possible since ak ^ A. Also, all ak must be positive roots 
because permutes the positive roots other than 7, for 7 G A. Then the loop in 1. must 
terminate because the height of ak decreases every step. 

Proposition 10 All Yk are elements ofU{n^) and we have an embedding M{fj,) ^ M(A) 
given by ^-^ Yr ■ vx. 

Proof. We write Si = For < A; < r we set = s^-fc ■ ■ ■ Si (so Wr = 1), and 
fJ^k = WkfJ', Xk = WkX. We claim that there is an embedding M{fik) ^ M{Xk) given by 

(— > Yk ■ v\^. First we look at the case /c = 0. Note that s.^ ■ ■ ■ Si{a) = /5 G A. Since for 
w in the Weyl group we have wSj3W~^ = 3^/3 we get Sa = si ■ ■ ■ SrSpSr ■ ■ ■ si = Wq^s^wq. 
Therefore /xq = woSa{X) = s^(Ao), and (Ao,/5^) = (A, si ■ ■ ■ Sr(/?)^) = (A,a^) = h. The 
case k = now follows by Lemma IHl 

Now suppose we have an embedding M{fj,k) ^ M{Xk) as above. Note that Wk+i = 
Sr-kWk and ak = Wka. Also fik = Wkfi = Xk — hak, and (A^, a^) = h, so that fik = Sa^{Xk)- 
We now apply Proposition |H1 (with u := /i^, r] := Xk, (3 := Pr-k)- We have Pr-k £ A and 
Pr-k 7^ ttfc as ak ^ A. Furthermore, m = {sr^k ■ ■ ■ Sifi, P^_k) = -(/i, Si ■ ■ ■ Sr-k-iiPr-kY) = 
ttr-k- In the same way n = {Xk, P^-k) — ~K-k- So by Proposition |H1 if we set 

Yk^. = y;Z:Y^yl-_:, 
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then we have an embedding M(/iA;+i) = M{sfi^_^iJ,k) ^ M{sp^_^Xk) = M{\k+i) by v^^^^ ^ 

Finally we note that = X, fir = fJ'- n 
It is possible to reformulate the algorithm in such a way that it looks more like the 
method from ^Hj- The construction described in JH] works as follows. Write Sa = 
■ ■ ■ Sa^^ , as a product of simple reflections. For 1 < A; < t set = {sa^^_^^ ■ ■ ■ ^a^X, a^). 
Then Y = ■ ■ ■y™'* is an element of U(n^) and we have an embedding M{fi) M(A) 
^ Y ■ vx. Now, using the same notation as in the description of the algorithm, the 
expression we get is 

ai ttr h br bi 

y/3i ■ ■ ■ yprtlpypr ■ ■ ■ y/3r 

As remarked in the proof of Proposition Sa = si - ■ ■ SrSpSr ■ ■ ■ si (where again we 
write Si = s^J. Furthermore, bk = (s^-i ■ ■ ■ SiA, /?^), h = {\,a^) = {sr ■ ■ ■ SiX, P^), 
= {sk+i ■ ■ ■ SrSisSr ■ ■ ■ SiX, So wc scc that our method is a special case of the con- 
struction in 18J. However, the difference is that we have an explicit method to rewrite the 
element above to an element of U{n~). By the next lemma the expression we use for Sa 
is the shortest possible (so we cannot do essentially better by taking a different reduced 
expression) . 

Lemma 11 The expression s^ = Si ■ ■ ■ s^S/js^ ■ ■ ■ Si obtained by the first step of the algo- 
rithm, is reduced. 

Proof. Set 7 = si{a) = a — m(3i, where m > 0. Then Sa = Ssi(-y) = Sis^si. By induction, 
the expression = S2 • • • SrSpSr ■ ■ ■ S2 is reduced. We show that l{sa) = ^(s^) + 2. For this 
we use the fact that the length of an element w of the Weyl group is equal to the number 
of positive roots that are mapped to negative roots by w. Write = A U {Pi}, where 
y4 = \ {Pi}. There is a positive root 5o ^ with s^Si(5o = Pi. Set S = {6 E A \ s^si6 < 
0} U {6o,Pi}. Then Sa maps all elements of S to negative roots. Since {'y,Pi) = —m < 0, 
also (/?i,7^) < 0, and hence s^(/3i) > 0. So all roots that are mapped to negative roots by 

are in A. Therefore, since si permutes A, there are £(7) roots 5 & A with s^si{5) < 0. 
We conclude that the cardinality of S is £(7) + 2. So £(«) > £(7) + 2, but that means that 
i{a) = £(7) +2. □ 

We can use the algorithm described in this section to construct general formulas for 
singular elements. More precisely, let 7 be a fixed root in the root system of g. Then by 
applying the formulas of Section |21 symbolically we can derive a formula that given arbitrary 
weights X, fi such that (A, 7^) G Z>o and = s^(A) produces an element Y G U{n^)\-^ 
such that Y ■ vx is singular. We illustrate this with an example. 

Example 12 Suppose that g is of type ^3. We use the same basis of as in Example 
El We consider the root a + /3 + 7. Let A = (ri, r2, r^) be such that h = ri + r2 + is a. 
positive integer. A reduced expression of Sa+i3+-y is SaSps^SisSa- The corresponding element 
of U{n-) is 
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First we have 

k=0 V / V / 

Now to obtain the formula for Y we have to apply Lemma El three times (and Lemma|21a 
few times), to obtain 

- 1 E g D-r^:) ^0 (7) ') (:) C 

hV}o]f]^A-l.-s h-k-s h-k-t s-t k-l l+t 

Table Q contains a few running times of the implementation of this algorithm in GAP4. 
The root 7 is in each case the highest root of the root system. The length of a formula is 



type 


length 


time (s) 


Ae 


29 


0.2 


De 


109 


1.6 




316 


2.9 




2866 


26.3 


Es 


> 10556 


00 



Table 1: Running times for the computation of a formula for a singular vector. 



the number of summations it contains (so the length of the above formula for ^3 is 4). The 
computation for Eg did not terminate in the available amount of memory (lOOM). When 
the program exceeded the memory, the expression contained 10556 summations. 
Remark. It is also possible to use this method to obtain formulas for a fixed type, but 
variable rank. However, for that a convenient Chevalley basis needs to be chosen. We refer 
to [IH], Section 5, for the formula for An. 

Remark. We have chosen Q as the ground field, because it is easy to work with. However, 
from the algorithm it is clear that instead we can choose any field F of characteristic zero 
and construct embeddings of Verma modules with highest weights from FP. 



4 Composition of embeddings 

In this section we consider the problem of obtaining an embedding M{v) ^ M(A), where 
z/ = So,Sj3{\) < sp{X) < A. The obvious way of doing this is to set = s^(A) and obtain the 
embeddings M{v) '—>■ M(/i), M{fi) M(A) and composing them. This amounts to multi- 
plying two elements of U{n~). This then corresponds to an expression for SaSjj, which is not 
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necessarily reduced. The question arises whether in this case it is possible to do better, i.e., 
to start with a reduced expression for s^S/j = Sa^^ ■ ■ -Sai^, set = {sa,^ ^ ■ ■ ■So.^A,a)^), 
and rewrite Y = y^^ ■ ■ -y™"* to an element of U{n~). The next example shows that this 
does not always work. 

Example 13 Let $ be of type F4, with simple roots ai, . . . , ^4 and Cartan matrix 

/ 2 -1 \ 
-12-2 
0-12-1' 
\ 0-12/ 

Let a = «! + 0^2 + 20:3 and /3 = ai + 2a2 + 2a3 + a4^. We abbreviate a weight aiAi + - ■ ■ + a4A4 
by (ai, a2, 03, 04). Set A = (|, — |, |, 0) and z/ = (— |, — i, — |, 2). Then z/ = SaS/3(A). Write 
Si = Sq,-. Then a reduced expression of SaSp is 

S1S2S1S3S2S1S3S2S4S3S2S1S3S2. 

We get 

121.53 45 ill _i _1 

y = Vi 1/2 Vi yi y2 yivi 2/2 Vivi y? % ' ?/2 ' • 

And I do not see any direct way to rewrite this as an element of U{n^). 

In general we have to obtain the embedding by composition. In this example set 
fi = Sf3{\) = A - /3 = (|, -|, |, 0). Then for the embedding M(/i) ^ M(A) we get 

2 4 2 1 _1 ill _i _1 

Yi = z/2 1/3 1/1 1/2 2/3 ""yiyiyiyiVs '2/2 

For the embedding M(z/) M{fi) we get 

y2 = 2/l2/3?/22/3 2/l • 

Then the product Y2Y1 will provide the embedding M{v) ^ M(A). 



5 Affine algebras 

In this section we comment on finding embeddings of Verma modules of affine Kac-Moody 
algebras. First we fix some notation and recall some facts. Our main reference for this is 

We let Q be the (untwisted) affine Lie algebra corresponding to g, i.e., 
with multiplication 

[f^m+aiK+bid, r®y+a2K+b2d] = y]+bint''®y-b2mt"'®x)+m6m,-nK.{x, y)K, 
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where m,n E Z, x,y & g, cti, 02, &2 £ Q and k( , ) is the Kilhng form on q. 

The Lie algebra q has a triangular decomposition g = © f) © n"*". Here n~ is spanned 
by the t™ © ?/j for m < 0, along with t"^ ^ Xi, and © hj for n < 0. The subalgebra f) is 
spanned by the © /ij and i^' and d. Furthermore, n"*" is spanned by the © Xi for m > 0, 
along with © ?/j and t" © hj for n > 0. 

The Verma module M(A) of highest weight A is defined in the same way as for g. As 
vector spaces M(A) = U{h'). Let a be a positive root of g. Then from [T2] we get that 
M(A — na) embeds in M(A) if and only if 2(A, a) = n{a, a), where tt, is a positive integer. 

Now if a is a real root with 2(A,a) = n{a,a), then we can construct a singular vector 
in U{h~)na by essentially the same method as in Section IHl The only difference is the 
algorithm for rewriting fl^~^af[, where r E Q, a E U{n^), and a basis element of 
n~. We need commutation relations y"^fj = fiV"^ + ■ ■ ■ 5 where y runs through the basis 
elements of n~. 

First of all, if fi = P © x^ for some a G and y = ® Xp for some /? G $ such that 
a + /3 G then set yma+np = © Xma+n/s- Set B = {yma+np \ma + n(3 E $}. Then 

B spans a subalgebra of g isomorphic to the subalgebra of g spanned by the corresponding 
Xma+np- The isomorphism is given by yma+nis ^ Xma+np- So we get the same formula as 
in the finite-dimensional case. 

Now suppose that = 0. Then j + A; < 0; so [fi, y] = P^^®ha, where ha = [x^, X-a\- 
In this case we use the following relation: 

r{P © XaY'\t''+' © K) - r(r - 1){P © XaY'^{t''+^' © 

which is easily proved by induction. If © X-a occurs with an exponent > 1 then we use 
this formula repeatedly. 
The last possibility is 

{t^ © hq){P © XaY = ® XaY{t^ ® K) + ® XaY'^{t^^^ ® Xa). 

Again, we use this formula repeatedly if © hq occurs with exponent > 1. 

Now we suppose that a = m6 is an imaginary root with (A, a) = (here 6 is the 
fundamental imaginary root). Then M(A — na) > M(A) for all positive integers n. In 
this case there are a lot of singular elements. One class of them is easily constructed. Let 
Ml, . . . , Mq, M^, . . . , m'^ be two basis of g, dual to each other with respect to the Killing form. 
For n > set 

q n 
i=l j=0 

Lemma 14 Suppose that (A, 6) = 0, then Sn-vx is a singular vector of weight n6 in M(A). 
Proof. From ^2), 12.8 we have the Sugawara operators 

Q 

meZ i=l 
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It is straightforward to see that Sn ■ vx = T^n ■ v\. Now K acts on M(A) as scalar 
multiphcation by — /i^, where is the dual Coxeter number. But also by |T2], Lemma 
12.8 we have for a; G g: 

[r ®x,T„„] = 2m(ir + /i^)(t'"-"®x). 

From this it follows that x ■ SnV\ = for < i < /, x G n^. Therefore 5*^ ■ f a is a singular 
vector. □ 

Lemma provides an infinite number of singular vectors. However, these are not the 
only ones. In J7] it is shown that for n > and 1 < i < I there are independent elements 
5** G U{n~) of weight nS, such that S"^ ■ f a is a singular vector. These 5"^ are constructed 
from the generators of the centre of U{g). In this construction the 5*^ correspond to the 
Casimir operator. However, with the exception of the Casimir operator, I do not know 
of efficient algorithms to construct the generators of the centre of U{q). For example, the 
explicit expressions given in [S] for a generator of the centre of degree s involve sums of 
(dimg)^ terms. So constructing generators of the centre of U{q) appears to be a very hard 
algorithmic problem in its own right. 

The conclusion is that we have efficient algorithms to construct an inclusion M(A — 
na) ^ M(A) if a is a real root, or when a is imaginary. However, in the last case there are 
many singular vectors for which at present we have no efficient way of constructing them. 



6 Constructing irreducible representations 

In pni; P- Littelmann proves a theorem describing a particular basis of the irreducible 
representations of g, using inclusions of Verma modules. Apart from giving a basis this 
result also allows one to construct the irreducible representations of q. In this section we 
first briefly indicate how this works, and then give some experimental data concerning this 
algorithm. 

The first ingredient of the construction is Littelmann's path method. Here we only give 
a very rough description of that method; for the details we refer to A path is a 

piecewise hnear function vr : [0, 1] MP, such that 7r(0) = 0. Such a path is given by two 
sequences Ji = (/xi, . . . , /i^) and a = (ao = 0, Oi, . . . , = 1), where the /ij G MP and the Oj 
are real numbers with = Oq < cti < • • • < = 1- The path vr corresponding to this data 
is given by 

s-l 

7i{t) = {t- a^_i)/is + - for a^-i <t <as. 

4 = 1 

Let A be a dominant weight. Then the path tta is given by the sequences (A) and (0, 1), 
i.e., it is the straight line from the origin to A. For a G A there is a path-operator fa. 
Given a path vr, /a(7r) is a new path, or 0. Set P(A) = {/a^^ ■ ■ ■ fai^ (tta) | ^ > 0, ai. G A}, 
and let ^(A) be the irreducible g-module with highest weight A. Then from PH] we 
have that the endpoints of the paths in P(A) are weights of V^(A) and the number of paths 
with endpoint is equal to the dimension of the weight space in V{X) with weight /i. 
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Let 71 G -B(A) be given by the sequences (/ii, . . . , fir) and (ao = 0, ai, . . . , = 1). Set 
Hr+i = A and z/j = Oj/Xj and ?7j = aifii^i for 1 < i < r. Then it can be shown that 
M^Vi) ^ M{rji). Let 0, e [/(n~),,._y- be such that • Vr,. is a singular vector. Then 
set 6^ = 6i ■ ■ - Qr- The element 0^ G t/(n~)A-7r(i) is determined upto a multiplicative 
constant. 

Now in PHI an inclusion B{mX) ^ B{n\) is described for m < n. With this inclusion 
we can view B[mX) 8bS clS 3i subset of B{nX). Furthermore, i?(A, cxo) denotes the union of 
all B{mX) for m > 1. Write A = riiXi + ■ ■ ■ + riiXi, and let /(A) be the left ideal of U{n~) 
generated by the elements y^'~^^ for 1 < i < I. Then V^(A) = M(A + p)/I{X) ■ v\, where 
p = Ai + ■ ■ ■ + A;. Now from (TH] we have the following result. 

Proposition 15 Suppose that all Ui > 0. Then the set {67^ | tt G -B(A, 00), tt ^ B{X)} is 
a basis o//(A). 

(If some rij = then there is a similar result, which we will omit here, cf. [T^.) 

In order to construct and work with the quotient M(A + p)//(A), we need a basis of 
/(A). If A — /i is not a weight of ^(A), then J(A) fl ?7(n")^ = ?7(n~)^. So we only need 
bases of the spaces /(A) fl f/(n~)^ where A — /x is a weight of V{X). By the above theorem 
we can compute those bases by first computing paths n G B{X, 00) with 7r(l) = A — /i, and 
then constructing the corresponding 0,^. We call this algorithm A. 

In Table El the running times are given of algorithm A on some sample inputs. Also 
listed are the running times of the algorithm described in :9|, which uses a Grobner basis 
method to compute bases of the spaces /(A) fl t/(n~)^. We call it algorithm B. In order 
to fairly compare both algorithms, the output in both cases consisted of the representing 
matrices of a Chevalley basis of g. 



type 


X 


dim\/(A) 


H inclusions 


time A (s) 


time B (s) 


A2 


(2,2) 


27 


64 


1.0 


1.3 


A2 


(3,4) 


90 


296 


2.3 


5.0 


A2 


(5,5) 


216 


788 


6.5 


14.2 


A; 


(1,1,1) 


64 


897 


16.6 


6.0 


A3 


(2,1,1) 


140 


2834 


56.4 


15.0 


A3 


(2,1,2) 


300 


7837 


178.4 


40.0 


B2 


(2,2) 


81 


807 


10 


6 


B2 


(3,3) 


256 


3330 


56 


23 


B2 


(4,4) 


625 


9502 


347 


79 



Table 2: Running times (in seconds) of the algorithms A and B for the construction of 
V{X). The fourth column displays the number of inclusions of Verma modules computed 
by algorithm A. The ordering of the fundamental weights is as in 

We see that for type ^2, algorithm A competes well with algorithm B. However, for 
the other types considered this is not the case. In these cases huge numbers of inclusions 
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of Verma modules have to be constructed, which slows the algorithm down considerably. I 
have also tried to construct ^(A) for A = (1, 1, 1), and g of type B^. But algorithm A did 
not complete this calculation within the available amount of memory (lOOM). 
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